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The combined effect of nonmagnetic and magnetic defects and impurities on critical temperatures of 
superconductors with different gap anisotropy is studied theoretically within the weak coupling limit of 
the BCS model. An expression is derived which relates the critical temperature to relaxation rates of 
charge carriers by nonmagnetic and magnetic scatterers, as well as to the coefficient of anisotropy of the 
superconducting order parameter on the Fermi surface. Particular cases of d-wave, (s + d)-wave, and 
anisotropic s-wave superconductors are briefly discussed. 

This paper is motivated by conflicting experimental results concerning the symmetry of the superconducting order 
parameter A(p) in high-temperature superconductors (HTSCs) and the suppression of the critical temperature Tc 
of HTSCs by defects and impurities. Indeed, while the majority (though not all) of experiments support the d-wave 
superconductivity in HTSCs , the observed degradation of Tc by impurities or radiation- induced defects |^ is more 
gradual than predicted theoretically for d-wave superconductors . 

To resolve this contradiction, a number of suggestions have been made, including anisotropic s-wave symmetry of 
A(p) 1^, momentum dependence of impurity scattering js), strong coupling effects resulting in crossover from Cooper 
pairs to local bosons 1^, etc. Note, however, that theoretical analysis of Tc degradation by defects and impurities 
is usually restricted to the speciflc case of spin-independent scattering potential, i.e., to the case of nonmagnetic 
scatterers only. Meanwhile a lot of experiments give evidence for the presence of magnetic scatterers (along with 
nonmagnetic ones) in non-stoichiometric HTSCs, e.g., in oxygen-deficient, doped or irradiated samples |^]. 

The goal of this paper is to work out a theoretical framework for a description of combined effect of nonmagnetic 
and magnetic scatterers on Tc of a superconductor with anisotropic A(p) (in what concerns an isotropic s-wave 
superconductor, its Tc is insensitive to nonmagnetic defects while the Tc suppression by magnetic defects is given 
by a well-known Abrikosov-Gor'kov theory 1^). We use the weak coupling limit of the BCS model for superconducting 
pairing and the Born approximation for impurity scattering. In what follows, we do not specify the microscopic 
mechanism of superconductivity. We set h = = ^ throughout the paper. 

The Hamiltonian of a superconductor containing both nonmagnetic and magnetic scatterers reads 



where ^(p) = e(p) — /i is the quasiparticle energy measured from the chemical potential, U{p, a; p', a') is the matrix 
element for electron scattering by randomly distributed impurities (defects) from the state (p', a') to the state (p, a), 
and V{p,p') is the BCS pair potential. 

We assume for simplicity that electron scattering is isotropic in the momentum space, the amplitude of the scattering 
by an isolated nonmagnetic (magnetic) scatterer being m„ (u,n). Then the relaxation times r„ and r™ are given by 
the standard "golden rule" formulas 



where c„ and Cm are the concentrations of scatterers, A^(0) is the density of electron states at the Fermi level. Note 
that the commonly accepted expression for is J^S'(S'-|-l)/4, where J is the energy of electron-impurity exchange 

interaction, S is the impurity spin. 

In order to account for anisotropy of the superconducting state, we assume a factorizable pairing interaction of the 
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where n = p/p is a unit vector along the momentum. Then the order parameter A(p) is 

A(p) = Ao0(n), 



(4) 



where Aq depends on the temperature. Thus the function 0(n) specifies the anisotropy of A(p) in the momentum 
space {4>{n) = 1 for isotropic pairing). The self-consistent equation for A(p) can be derived by means of Green's 
functions technique (see, e.g., 0]). It is as follows: 



A(p) = - ^F(p,p')(a_p,xVT) = -TY,Y.^{p,p') 
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where uj = 7rT(2n + 1) are Matsubara frequencies, and the equations for Ai^{p) and ui' are 
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Uj' = uj - i(c„|u„p + c„|u„p) ^ ■ 



tLo'+ap') 



-C.'2+^2(p,)+|A^(p')|2' 

Since A(p) = at T = Tc, in the case T ^ Tc we have from (|^), (Q), taking into account: 
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Lli' =Lli + ^ (l/r„ + l/r„) sign{uj), 
where angular brackets (...) stand for the average over the Fermi surface (FS): 
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Substituting (^) and (^) in (^ and taking (^ into account, we have after rather simple but time consuming algebraic 
transformations: 



\uj\ + \ (l/r„ + l/r„) 
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Here Tco is the critical temperature in the absence of impurities and defects (at l/r„ — 1/t„i — 0). At this stage it is 
convenient to introduce the coefficient x of anisotropy of the order parameter on the FS |10|, B 
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For isotropic s-wave pairing we have A(p) = const on the FS; therefore, (A(p))2 = (A2(p)), and x = 0. For a 
superconductor with d-wave pairing we have x = 1 since (A(p)) — 0. The range < x < 1 corresponds to anisotropic 
s-wave pairing or to mixed {d + s)-wave pairing. The higher is the anisotropy of A(p) (e.g., the greater is the partial 
weight of a d-wave in the case of mixed pairing), the closer to unity is the value of x- 
Making use of the definition ([l2|) and the formula pi| 
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where ^I' is the digamma function, we obtain from (|11| 
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In two particular cases of (i) both nonmagnetic and magnetic scattering in an isotropic s-wave superconductor (x = 0) 
and (a) nonmagnetic scattering only in a superconductor with arbitrary anisotropy of A(p) (1/t„j = 0, < x 1)j 
the Eq.(nJ) reduces to well-known expressions M, mn 
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respectively. 

The Eq. (|lj) is obviously more general than Eqs.(15) and (16) which are commonly used for the analysis of 
experimental data on Tc suppression by defects and impurities in HTSCs fl^ . In fact, making use of Eq. ([l5| ) or 
Eq.(|l6|) one assumes a priori that either (i) the order parameter in HTSCs is isotropic in the momentum space or 
(ii) the magnetic scatterers in HTSCs are completely absent. The latter assumption is often supplemented with a 
speculation about pure d-wave symmetry of A(p) fl^ (i.e., one intentionally restricts himself to the case x = 1 instead 
of attempts to extract the value of x from the experiment). In our opinion, the experimental dependencies of Tc versus 
impurity (defect) concentration or radiation dose should be analyzed within the framework of the theory presented 
above, see Eq.(|lJ). One should not guess as to the degree of anisotropy of A(p) and the type of scatterers, but try to 
determine the value of x f^nd relative weights of magnetic and nonmagnetic components in electron scattering through 
comparison of theoretical predictions with available or specially performed experiments. 

Now let us consider the limiting cases of weak and strong scattering {Tco — Tc << TcO and Tc — > respectively). 
At 1 / AnTcoTn < < 1 and 1 / AirTcoTm < < 1 (weak scattering) one has from ( p^ ) : 
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In particular cases [i] and {ii) considered above, Eq.(|T^ reduces to well-known expressions 
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for initial Tc reduction by magnetic (at X — ^) or nonmagnetic (at arbitrary value of x) scatterers respectively. 

In what concerns the strong scattering limit, we recall that in the BCS theory, nonmagnetic scattering alone is 
insufficient for the not-d-wave superconductivity (0 < x < 1) to be destroyed completely |l^; at 1/t,„ = 0, the value 
of Tc asymptotically goes to zero as 1/t„ increases (whereas Tc of a d-wave superconductor with x — ^ vanishes at a 
critical value l/r,"; — ttTcoM ~ l-764rco, where 7 = e*-^ « 1.781, C is the Eiler constant). On the other hand, magnetic 
scattering in the absence of nonmagnetic one (I/t^ = 0) is known to suppress the isotropic s-wave superconductivity 
with X = at a critical value 1/r^ = iiTco/^-f k. 0.882rco [§. 

Based on the Eq.([T^), it is straightforward to derive the general condition for impurity (defect) suppression of 
Tc of a superconductor having an arbitrary anisotropy coefficient x and containing both nonmagnetic and magnetic 
scatterers: 
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where r^^ j is the critical value of the effective relaxation time Tcf / defined as 



Teff 



(20) 
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From Eqs. (gg) and (|2l|) one can see that increases monotonically with both l/r„ and l/r^ at any value of 

X, with the exception of the case x = when I/te// doesn't depend on l/r„, see (21). If x is close to unity (strongly 
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anisotropic A(p)), then ^/T^ff ~ 1/t„ + l/rm, i.e., the contribution of nonmagnetic and magnetic scattering to 
pair breaking is about the same. If x « 1 (almost isotropic A(p)), then '^/reff w I/t™, i.e., Te// is determined 
primarily by magnetic scattering. The higher is the anisotropy coefficient x, the greater is the relative contribution 
of nonmagnetic scatterers to degradation with respect to magnetic ones. If nonmagnetic scattering is absent 
(1/t„ = 0), then ^/Teff = 1/tv„ at any value of x- 

We note however that while the concept of the effective relaxation time Te// can be used for evaluation of the 
critical level of nonmagnetic and magnetic disorder, it is not possible to express Tc in terms of Tg// in the whole range 
Q < Tc < Tco, see (p^). In other words, the combined effect of nonmagnetic and magnetic scattering on Tc cannot 
be described by a single universal parameter depending on the values of r„, t^, and x- For example, I/te// = at 
I/tto = and < % < 1 no matter what the value of 1/t„ is. On the one hand, as follows from (|20|), the zero value 
of in this case points to the fact that in a BCS superconductor with not-d-wave symmetry of A(p) the critical 

level of disorder cannot be reached in the absence of magnetic scattering, in accordance with . On the other hand, 
the zero value of 1 /re/ / obviously doesn't imply that Tc of a not-d-wave superconductor is completely insensitive to 
nonmagnetic scatterers at l/r™ = and < x < li see (|4|). Hence, while the quantity I/t^j^ characterizes the 
critical strength of impurity (defect) scattering corresponding to Tc — 0, the quantity XjTcff (when it is less than 
l/r^^j) doesn't determine the value of Tc unequivocally. 

Based on Eqs. (|20| ) and (|2ll), it is possible to derive the following expression for the critical value of 1/t„ in the 
presence of magnetic scattering: 



This expression is valid as long as l/r™ < 7r2^^^Tco/7 since otherwise the superconductivity is completely suppressed 
solely by magnetic impurities. The value of decreases as 1/t,„ increases at constant x or as x increases at 
constant l/rm,. The finite value of in the presence of magnetic scatterers could reconcile the experimentally 
observed disorder-induced suppression of Tc of HTSCs below 4.2K [| with theories of not purely d-wave symmetry 
of A(p) in HTSCs, e.g., anisotropic s-wave symmetry or mixed (d + s)-wave symmetry. 

In conclusion, the results obtained provide the basis for evaluation of the degree of anisotropy of the superconducting 
order parameter (and hence its possible symmetry) as well as the type of scatterers (magnetic or nonmagnetic) in high- 
Tc superconductors through careful comparison of theoretical predictions with the experiments on impurity-induced 
and radiation-induced reduction of the critical temperature. We hope that the present paper will serve as a stimulus 
for experiments on combined effect of nonmagnetic and magnetic scattering in the copper-oxide superconductors. 

This work was supported by the Russian Foundation for Fundamental Research under Grant No 97-02-16187. The 
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Note added in 'proof. After submission of this paper I became aware of a similar study by A.A.Golubov and 
I.I.Mazin [Phys. Rev. B 55, 15146 (1997)] which generalizes Abrikosov-Gor'kov solution to the case of a multiband 
superconductor with interband order parameter anisotropy. 
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